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[57] ABSTRACT 

The columnar electron content of the ionosphere be- 
tween a spacecraft (10) and a receiver (12, 14) is mea- 
sured in realtime by cross correlating two coherently 
modulated signals transmitted at different frequencies 
(Li, L 2 ) from the spacecraft to the receiver using a cross 
correlator (16). The time difference of arrival of the 
modulated signals is proportional to electron content of 
the ionosphere. A variable delay (21) is adjusted relative 
to a fixed delay (20) in the respective channels (Li, L 2 ) 
to produce a maximum at the cross correlator output. 
The difference in delay required to produce this maxi- 
mum is a measure of the columnar electron content of 
the ionosphere. A plurality of monitoring stations (31, 
32, 33) and spacecrafts (GPS satellites 34, 35, 36) are 
employed to locate any terrestrial event that produces 
an ionospheric disturbance. 

7 Claims, 4 Drawing Figures 
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METHOD AND APPARATUS FOR CALIBRATING 
THE IONOSPHERE AND APPLICATION TO 

SURVEILLANCE OF GEOPHYSICAL EVENTS 

ORIGIN OF INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA contract and is sub- 
ject to the provisions of Section 305 of the National 
Aeronautics and Space Act of 1958, Public Law 85-568 
(72 Stat. 435; 42 USC 2457). 

BACKGROUND OF THE INVENTION 

This invention relates to a technique for calibrating 
the ionosphere in respect to columnar electron content 
on a realtime universally applicable basis, and to appli- 
cation of that calibrating technique to surveillance of 
large-scale geophysical events which modify the iono- 
sphere. 

The ionosphere has been extensively explored in 
nonrealtime. The earliest technique involved the 
Ionosonde utilizing a pulsed transmitter that sends radio 
signals vertically upward while slowly sweeping the 
radio frequency. 

Pulses reflected by the ionosphere are received and 
recorded; the elapsed time between pulsed transmission 
and reception can be interpreted as an apparent distance 
to the point of reflection. This sounds only the bottom 
side of the ionosphere. The resulting presentation is a 
curve of apparent height vs. frequency. This technique 
is still profitably and extensively used today. In recent 
years, many rockets have flown through the ionosphere 
to record data, and satellites have orbited above the 
ionosphere for similar purposes. An ionosonde in a 
satellite has been utilized to explore the topside of the 
ionosphere (that portion above the region of maximum 
electron concentration). 

The Ionosonde technique applied at ground stations 
showed that the ionosphere is structured. It was first 
thought that discrete layers were involved, commonly 
referred to as D, E, Fi, F 2 . However, rocket measure- 
ments have shown that the “layers” merge with one 
another to such an extent that they are now normally 
referred to as regions. A vertically incident radio wave 
is reflected from the ionosphere at that level where 
natural frequency of the plasma equals the radio fre- 
quency. This procedure makes it possible to identify the 
electron concentration at the point of reflection in terms 
of the radio frequency according to the relation 

n c =1.24Xl0 10 / (1) 

where n e is the electron concentration per cubic meter 
and f is the radio frequency in megahertz. 

The D region is the lowest ionospheric region ex- 
tending approximately from 60 to 85 kilometers. The D 
region is susceptible to disturbance particularly due to 
x-ray bursts. The sudden ionospheric disturbances 
(SID’s) occur over the daylight hemisphere with some 
solar flares; it can be shown that SID’s also occur be- 
cause of acoustic pressure, atomic blasts (including un- 
derground), volcanic eruptions, earthquakes which 
cause acoustic percussion, and large tropospheric 
storms. 

VHF Faraday rotation has been used to estimate the 
charged particle count of the ionosphere, but this Fara- 
day effect depends upon Earth’s magnetic field and no 
measurable effect exists for electrons out of the mag- 
netic field influence. Therefore, this technique cannot 


inherently measure the total electron content without 
modeling associated parameters. 

A surprisingly large number of applications for accu- 
rate ionospheric measurements exist. The ionosphere 
5 represents a significant error source for high precision 
radio metric tracking data. For example, 1 X 10 17 el/m 2 
zenith columnar content during night hours results in 12 
cm of baseline length increase on an Astronomical 
Radio Interferometric Earth Surveying (ARIES) geo- 
10 detic measurement at S-band on a 336 km distance. 
Spacecraft tracking at S-band and to a lesser degree at 
X-band by a factor of 13, also suffers from errors caused 
by free electrons in the radio path which are substan- 
|5 tially the Earth’s ionosphere. The time variable iono- 
spheric electron columnar content in the tracking sta- 
tion to spacecraft line of sight causes variations which 
are highly correlated with the Doppler effects caused 
by changes in spacecraft right ascension and declination 
2 0 positions. 

Less prosaic applications also exist because of the 
effect of troposphere to ionosphere dynamic coupling. 
Disturbances in the troposphere propagate into the 
ionosphere and can cause propagating ionospheric 
25 waves. Speculations exist concerning impulsive input in 
the troposphere which has manifestation in the iono- 
sphere. For example, a tsunami event in the ocean 
causes an acoustic impusles which propagates into the 
ionosphere. The March 1980 Mount St. Helens eruption 
30 was very clearly evident in VHF Faraday rotation 
ionospheric monitoring at Goldstone, Calif. That Fara- 
day rotation instrumentation has a sensitivity of approx- 
imately 3 X 10 15 el/m 2 on 1 minute samples and an ambi- 
guity of 4 X 10 16 el/m 2 . Other impulsive events of inter- 
35 est are nuclear explosions, especially underground ex- 
plosions in which seismic hiding methods might be 
employed to confuse seismic detection of such events. If 
conducted underground, the explosive impulse could be 
transmitted to the surface and result in a tropospheric 
40 impulse which may be detectable as an ionospheric 
columnar electron content modification. 

More direct influences in the ionosphere will also be 
detectable such as particle interactions, principally solar 
induced, but perhaps the result of ionospheric heating 
45 by input electromagnetic energy or disturbances caused 
by hot gases expelled by missiles during their boost 
phase or perhaps the testing of particle beam weapons. 

Large scale weather systems in the troposphere and 
stratosphere probably also effect ionospheric patterns 
providing yet another opportunity to study global mete- 
orological processes. Combining these ionospheric ob- 
servations with the direct range and Doppler measure- 
ments from high precision geodetic receivers will yield 
55 small scale variations along the line of sight which are 
caused by tropospheric index of refraction changes. 
The water vapor and dry tropospheric components are 
responsible for these index changes which can be de- 
tected with millimeter path equivalent sensitivity once 
60 ionospheric effects have been removed. 

It is thus evident that if the ionosphere is continually 
calibrated, any sudden change in the ionosphere will 
signal a terrestrial event within seconds. It is important 
to detect some events in realtime, such as a tsunami 
65 event in order to alert those in the path of these destruc- 
tive sea waves, or the launching of an intercontinental 
missile in order to launch an intercepting missile. Other 
events, such as an underground atomic explosion of an 
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experimental nature may not be as critical, but the abil- 
ity to detect it in realtime would still be important. 

The problem then is one of providing a technique for 
continual calibration of the ionosphere in realtime on at 
least a regional scale, but preferably on a global scale. 
This capability has not heretofore been available. Since 
the localized ionosphere is modified by terrestrial 
events because of troposphere acoustic coupling and 
direct atomic interaction, such as is due to hot gases 
expelled from missiles in their boost phase or plasma 
clouds in front of ballistic reentry vehicles or particle- 
beam impingng on the ionosphere, a plurality of iono- 
sphere surveillance stations could then be clustered 
about any strategic area to provide passive continuous 
monitoring of that area for military applications as well 
as civil, commercial and geodetic applications while 
being compatible with existing active and passive local 
surveillance systems. 

SUMMARY OF THE INVENTION 

In accordance with the invention, the columnar elec- 
tron content of the ionosphere between a spacecraft and 
a receiver at a ground station is measured in realtime by 
cross correlating two coherently modulated signals 
transmitted at different frequencies from the spacecraft 
to the receiver. Both frequencies are coherently modu- 
lated by noise, or pseudo random noise (PRN) by the 
transmitter. A delay in the transmission of one fre- 
quency signal relative to the other is a measure of the 
integrated ion concentration of the ionosphere. Any 
large and sudden change in the ion concentration could 
indicate a large-scale terrestrial event. By monitoring 
any changes in the delay in transmission of one fre- 
quency signal versus the other from a plurality of space- 
craft relative to a plurality of receivers, the approximate 
location of a terrestrial event causing a perturbed region 
of the ionosphere can be determined from the raypath 
directions of the receivers from the spacecraft transmit- 
ters. The present invention thus provides a means for 
continuously monitoring the ionosphere at a plurality of 
stations which continuously measure the columnar elec- 
tron content of the ionosphere such that the location of 
any sudden ionospheric disturbance occurring between 
stations can be determined in realtime, constituting a 
global surveillance system. 

The novel features that are considered characteristic 
of this invention are set forth with particularity in the 
appended claims. The invention will best be understood 
from the following description when read in connection 
with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic diagram illustrating the tech- 
nique for calibrating the ionosphere in accordance with 
this invention. 

FIG. 2 illustrates in waveforms A, B and C a PRN 
code for modulation on a carrier, the carrier modulated 
by the PRN code, and the PRN code detected from the 
carrier, respectively. 

FIG. 3 illustrates an exemplary cross correlator for 
the technique of FIG. 1. 

FIG. 4 illustrates an application of the technique for 
calibrating the ionosphere surveillance of large-scale 
geophysical events by passive realtime monitoring of 
the charged particle content of the ionosphere. 


DESCRIPTION OF PREFERRED 
EMBODIMENTS 

Referring to FIG. 1 of the drawings, the ionosphere 
5 calibration system exploits the fact that a Global Posi- 
tioning System (GPS) satellite 10 simultaneously broad- 
casts two L-band channels, Li at 1.575 GHz and L 2 at 
1.227 GHz, with a 20 MHz bandwidth pseudorandum 
noise (PRN) coherently modulated on both channels. 
10 These two PRN modulations are emitted in phase coin- 
cidence, but the L 2 channel modulation will be delayed 
by the ionosphere a greater amount than the 1.575 GHz 
modulation according to a frequency squared relation- 
ship. The time difference of arrival of the L 2 modulation 
15 relative to the Li modulation is directly proportional to 
the total columnar electron content along the observer- 
to-satellite path. It is precisely this dual Li and L 2 dis- 
persion relationship which was explicitly designed into 
the GPS for range calibration from its inception. How- 
20 ever, GPS user calibrations depend upon having a de- 
coder for the received Li and L 2 PRN sequences in 
order to correct the Li pseudo range measurement for 
typical GPS applications. Where only an ionospheric 
measurement is desired, it is unnecessary to explicitly 
25 know the PRN sequence. All that is needed is exploita- 
tion of the fact that the same sequence arrives on the L 2 
channel as arrived first on the Li channel, and to use 
cross correlation of Li and L 2 to determine this time 
difference. In this manner, full use is made of the 20 
30 MHz bandwidth transmissions without knowledge of 
the details of the PRN sequence. Information about the 
PRN codes is generally unavailable and the GPS codes 
may even be encrypted. 

An antenna 14 directed at the GPS is a conventional 
35 dual frequency antenna feeding two channels Li and L 2 
in the receiver, where channel Li is tuned to the fre- 
quency 1.575 GHz, and channel L 2 is tuned to the fre- 
quency 1.227 GHz. FIG. 2 illustrates in waveform A a 
typical PRN code used to modulate the channels Li and 
40 L 2 , and waveform B illustrates the nonreturn-to-zero 
modulation of the PRN code on one of the channels. 
The receiver amplifies and separates the modulated Li 
and L 2 signals into two channels which are then con- 
verted to digital form for cross correlation (or simply 
45 cross correlated in analog form) to determine the delay 
which must be introduced at the receiver to bring the 
two signals into coincidence, at which time the output 
of the cross correlator is maximized, or minimized, 
depending on the implementation of the cross correla- 
50 tor. 

An analog implementation of the cross correlator is 
illustrated in FIG. 3. Briefly, an amplifier 18 first ampli- 
fies the signals which are then separated into two chan- 
nels Li and L 2 , such as by filters at the inputs of the 
55 channels. The signals are converted to baseband in the 
channels, which is to say the modulations on the sepa- 
rate signals at frequencies Li and L 2 are obtained sepa- 
rately in each channel. The output of the channel L 2 is 
applied to an RF mixer 19 through a fixed delay 20 
60 while the output of the channel Li is applied to the RF 
mixer through a variable delay 21. The latter is adjusted 
for a maximum voltage output from the mixer, as indi- 
cated by an RMS voltmeter 22. The difference between 
the values of the fixed delay and the variable delay is a 
65 measure of the columnar electron content between the 
transmitter and the receiver. 

The variable delay may be implemented with a plu- 
rality of separate delay lines, each a known increment 
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longer than the other. A switch is then operated to 
connect each successive one in sequence until a maxi- 
mum output is produced from the mixer. Programmed 
digital control of the switch can provide for fast and 
automatic measurements to monitor the ionosphere 5 
continuously without the need for an operator. Alterna- 
tively, both the measurement circuits and the control 
can be implemented with digital circuits. 

Preliminary accuracy estimates of this satellite iono- 
spheric calibration technique for L-band are 10 16 elec- 10 
trons per meter squared absolute and perhaps 10 15 
el/m 2 for differential total electron contents which is 
the important quantity in radio interferometry data 
types. For S-band quasar interferometry applications, 
the station to satellite raypaths can be calibrated to an 15 
accuracy of approximately one centimeter with some 
accuracy degradation encountered in mapping to the 
quasar directions. At X-band, the calibration accuracy 
should be path equivalent to the submillimeter level. 
For applications to GPS geodetic missions, the iono- 
sphere should be calibrated to ah accuracy of 5 cm, of 
better, limited probably by the signal to noise conditions 
in the cross correlation processing of the 20 MHz band- 
width signals in the Li and L 2 channels. 

Conceptually, this same approach can be applied to 
the tracking of S/X-band downlink equipped spacecraft 
which require an uplink ranging code. The spacecraft 
ranging channel is only a frequency translating micro- 
wave repeater. Thus, whatever noise is present within 
the ranging channel bandwidth (i.e., 1.5 MHz baseband) 
is coherently modulated onto the S and X channels. By 
tracking the S and X channels, a measurement of the 
columnar electron content is made without the neces- 
sity for ranging in either channel. Several economies are 
possible with this innovation: 

1. Fewer kilowatts of power need be consumed in 
generating high power S-band Earth to spacecraft 
transmissions; 
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PRN sequence; it is sufficient to cross correlate the 
outputs of the two receiver channels. 

An application of this ionosphere calibration tech- 
nique will now be described with reference to FIG. 3. 
This application is for global monitoring of events such 
as missile launchings, atomic explosions, earthquakes, 
volcanoes, sun-spot activity, and other phenomena 
which affect the ionosphere or are influenced by the 
ionosphere. No system presently exists to accurately 
monitor the ionosphere continuously in realtime espe- 
cially without detailed GPS code knowledge, though 
this is potentially a powerful source for early warning 
and forecasting information for civil, commercial and 
military purposes. The GPS satellites, already globally 
available, can be utilized to furnish dual-frequency bea- 
con signals on a broad spatial basis. Other satellites to be 
added by the mid 1980’s can be utilized for higher spa- 
tial resolution. When fully implemented, the GPS cov- 
erage of the Earth will be complete such that at any 
location at any time, at least four satellites will be above 
the local horizon. Thus, four raypaths can be sampled at 
any arbitrary time. 

FIG. 4 is a schematic diagram illustrating the broad 
two-dimensional coverage of three ground stations 31, 
32 and 33 of an ionosphere radio surveillance system 
based on a plurality of two-frequency signal beacon 
sources, such as the L-band GPS satellites 34, 35 and 36. 
The two-frequency signals are received at the ground 
stations by dual-channel correlation receiver systems 
for columnar sampling and monitoring of the iono- 
sphere charged particle content, as described with ref- 
erence to FIGS. 1 and 2 for one station and one GPS 
satellite with reference to FIGS. 1 and 2. 

Ground stations 31, 32 ad 33 outlying a strategic area 
A normally monitor the dual-frequency signals emanat- 
ing from satellites 34, 35 and 36, establishing an average 
charged particle content along the raypaths through the 
ionosphere from the satellites to the antennas. A pertur- 


2. Longer microwave transmitter tube life in Earth bation of the charged particle content of the ionosphere 
stations through less high power use is possible; 40 caused by a volcanic eruption or even a missile launch- 

3. Smaller amounts of data need be processed com- ing, for example, can then be located in two dimensions 

pared to obtaining ranging at S-band to be differ- from the raypath directions of the antennas. Further- 

enced against ranging at X-band to reveal total more, any phenomenon traversing two or more 

electron content in the Earth to spacecraft path. raypaths; i.e., showing a similar but sequential perturba- 

In addition, during times of solar superior conjunction 45 tion of the ionosphere, can yield velocity as well as 
tracking of spacecraft, the S-band uplink becomes so direction of travel of the phenomenon, 
disturbed that the spacecraft receiver loses phase lock It will be appreciated that the ionospheric radio sur- 
and, as a result, can no longer coherently transpond the veillance system is a completely passive system which 

ranging code. Use of cross correlation charged particle cannot be detected and which is compatible with all 

calibration continues to function in this environment 50 other types of surveillance systems since it depends only 
because the spacecraft received S-band ranging channel on the reception of signals and interpretation of charged 

is only noise, but coherently modulated onto the S/X particle measurements through the ionosphere. Al- 

downlink, and is fully useable in a cross correlation though the preferred embodiment uses GPS satellites 

implementation to its full 1.5 MHz bandwidth. Conse- modulated by a PRN code, it is unnecessary for any 

quently, it should be understood that the present inven- 55 ground station to explicitly know the PRN code. All 
tion is not limited to GPS satellites, although the GPS that is needed is exploitation of the fact that the same 

satellites operating as currently planned with PRN sequence and/or digital transition rate arrives on the L 2 

modulation are well suited. channel as first arrived on the L; channel, and to use 


Once all GPS satellites are operational, there will cross correlation of Li and L 2 to determine this time 
always be at least four in view from any Earth location, 60 differene. In this manner, full use is made of the 20 MHz 
each coherently modulated on two L-band carriers. As bandwidth transmissions without knowledge of the 
these two carriers are received from GPS satellites at PRN sequence. 

different ground stations, each using a directional an- Although particular embodiments of the invention 
tenna, the difference in time of arrival of the two chan- have been described and illustrated herein, it is recog- 
nels at each antenna is a function of the delay due to the 65 nized that modifications and variations may readily 
ionosphere, i.e., the electron content of the column of occur to those skilled in the art. Consequently, it is 
space between the satellite and the ground station. As intended that the claims be interpreted to cover such 
noted hereinbefore, it is not necessary to decode the modifications and equivalents. 
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What is claimed is: 

1. A method of calibrating the electron columnar 
content in the space between a transmitter on a space- 
craft and a receiver on Earth comprising the steps of: 
transmitting from said transmitter two broadband sig- 
nals separated in frequency, each coherently modulated 
by pseudorandom noise; receiving in a direct path from 
said transmitter both of said signals at a receiver having 
two channels for obtaining said modulations on said two io 
signals; and cross correlating the modulation of one 
received signal to the modulation of the other received 
signal to determine the phase difference between the 
modulation of said one received signal to the modula- 
tion of said other received signal as a measure of elec- ^ 
tron content in the space between the transmitter and 
receiver. 

2. A method as defined in claim 1 wherein said cali- 
bration is carried on continuously in real time to moni- 20 
tor the electron columnar content for any change 
therein. 

3. A method as defined in claim 2 wherein said space- 
craft is at a known position relative to the Earth and 
said receiver is at a known position on Earth and 
changes in the electron columnar content between the 
Earth and said spacecraft is indicative of a terrestrial 
event of significance, including the step of detecting any 
sudden change in electron columnar content in the 
known path between said transmitter and said receiver. 

4. A method as defined in claim 3 used a plurality of 
receivers and transmitters, each transmitter on a differ- 
ent spacecraft of known positive relative to the Earth, 
and each receiver being in a known position on the 35 
Earth for monitoring the columnar electron content 
between receivers and transmitters, each receiver oper- 
ating the same as all other receivers with signals re- 
ceived from a transmitter, including the step of locating 
the source of sudden change in the electron content in 
the space between said plurality of transmitters and 
receivers in two dimensions from the known paths of 
the receivers to the transmitters. 

5. A method as defined in claim 4 wherein a plurality 45 
of spacecraft transmitters are in a position to be in line 
of sight communication with each receiver at all times. 
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6. Apparatus for calibrating the electron columnar 
content between a transmitter in space and a receiver on 
Earth comprising 

means for transmitting from said transmitter two 
broadband signals separated in frequency, each 
coherently modulated with pseudorandom noise, 
means at said receiver for receiving said transmitted 
signals along a direct path from said transmitter, 
means for cross correlating said modulation of one 
received signal with the modulation of the other 
received signal at said receiver, including 
means for separating said two received signals into 
two channels converted to baseband, 
means for adjustably delaying the baseband signal 
in one channel relative to the other, 
means for cross correlating said delayed baseband 
signals, and 

means for determining the delay of said first chan- 
nel signal producing the maximum output of said 
cross correlating means as a measure of the co- 
lumnar electron content between said transmit- 
ting means and said receiving means. 

7. A system for detecting the occurrence of a terres- 
trial event that disturb the ionosphere by monitoring 

25 electron content of the ionosphere using a plurality of 
stations and spacecrafts at known locations, each station 
operating the same as all other stations with signals 
received from a spacecraft, comprising 
means at each spacecraft for transmitting two signals 
at different frequencies coherently modulated with 
pseudorandom noise, 

means at each station for receiving said two signals in 
a direct path from a spacecraft, 
means for converting each of said two signals re- 
ceived to baseband, 

means for adjustably delaying one baseband signal 
relative to the other, 

means for cross correlating said two baseband signals 
as said one is delayed relative to the other to deter- 
mine the delay of said one relative to the other 
which produces a maximum cross correlation out- 
put as a measure of said electron content, and 
means for locating the position of the source causing 
disturbance of the electron content of the iono- 
sphere from the known position of stations and 
spacecrafts detecting a disturbance. 

***** 
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